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Abstract 

We introduce the concept of kinetic equations representing a natural extension 
of the more conventional notion of a kinetic relation. Algebraic kinetic relations, 
widely used to model dynamics of dislocations, cracks and phase boundaries, link 
the instantaneous value of the velocity of a defect with an instantaneous value of the 
driving force. The new approach generalizes kinetic relations by implying a relation 
between the velocity and the driving force which is nonlocal in time. To make 
this relations explicit one needs to integrate the system of kinetic equations. We 
illustrate the difference between kinetic relation and kinetic equations by working 
out in full detail a prototypical model of an overdamped defect in a one-dimensional 
discrete lattice. We show that the minimal nonlocal kinetic description containing 
now an internal time scale is furnished by a system of two ordinary differential 
equations coupling the spatial location of defect with another internal parameter 
that describes configuration of the core region. 

Keywords: martensitic phase transitions, lattice dynamics, kinetic relations, kinks, de- 
fects, quasicontinuum models, dispersion, nonlinear waves 

1 Introduction 

Kinetic relations attributing a particular value of velocity to a given value of the driving 
force are widely used in continuum mechanics as a constitutive description of such lattice 
defects as phase boundaries, dislocations and cracks (see the reviews [21 [IHl [26l [28] ) . These 
algebraic relations form independent postulates that serve as closing conditions specifying 
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singular solutions in the classical continuum theories. Since kinetic relations replace the 
detailed modeling of the core regions of the defects, they represent a condensed description 
of the complex physical behavior at the microscale. In practice, kinetic relations are either 
taken from experiment or deduced from the solutions of auxiliary microscale problems. 
Characteristically, these auxiliary problems always assume constant values of both the 
macroscopic velocity of the defect and of the corresponding driving force. 

Due to the implicit assumption that the defect is in a steady motion, the kinetic 
relations-based description misses the details of a nonsteady internal dynamics of the core 
region. In particular, the internal pulsations originating from the defect interaction with 
localized micro-inhomogeneities become averaged out. This is in contradiction with the 
presence of accelerated motions of defects at all scales, as revealed, for instance, by the 
power law acoustic emission accompanying plasticity, martensitic phase transitions and 
fracture (e.g. [21]). 

To partially recover the missing information we propose in this paper to replace alge- 
braic kinetic relations by differential kinetic equations. The aim of these equations is to 
capture the transient phases of the defect evolution in response to nonsteady driving. In 
the language of constitutive theory we propose to replace the instantaneous rheological 
relations on the phase boundary by nonlocal memory functionals originating from a local 
description in terms of internal variables. Such extension brings into the conventional 
theory an internal time scale and allows one to deal with the so-called "rate effects" . An 
example of a similar development is provided by the rate-dependent constitutive laws in 
the theory of friction, where the set of internal state variables is also assumed to satisfy 
differential constitutive relations (e.g. [311 [32]). 

In order to show the possibility of a systematic derivation of kinetic equations from 
a micromodel we consider a prototypical overdamped defect moving in a lattice. We de- 
velop a low-parametric description of the internal dynamics of this defect involving some 
specially selected internal variables which characterize the structure of the core region. 
Our approach can be viewed as an example of a quasicontinuum method in the sense of 
[H EZl ES] whose goal is to match the macroscopic continuum description outside singular- 
ities with a more detailed atomistic resolution of the core regions. The zero-order model 
of this type produces kinetic relations and implies instantaneous and universal response 
of the core region to external perturbations (autonomous core region). The first-order ap- 
proximation leads to kinetic equations which already capture some nonuniversal features 
of the microscopic dynamics. For consistency of the two approximations, kinetic equa- 
tions must of course reduce to the kinetic relation when the variation of the macroscopic 
parameters is sufficiently slow. 

To illustrate the main idea of our approach consider a toy model describing an over- 
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damped dynamics of a configurational point in a one-dimensional energy landscape: 



dv 



du 



(1) 



Here z/ is a variable defining the microstate of the system and G{t) is a sufficiently slow 
varying macroscopic driving force which depends on slow time t. The fast time is defined 
as r = t/e, where e << 1 is a small parameter. We assume that the gradient of the 
energy landscape is periodic in u. The method of kinetic relations postulates 

the existence of an algebraic relation between the driving force G and the macroscopic 
velocity 



Such relation, which we write as V{t) = r{G{t)), can be computed explicitly if G is a tilt 
of the energy landscape 



In the realistic situations the microscopic description of the type ([H) is too complex because 
it involves a huge number of variables. In contrast, the macroscopic description (jl]) is 
too schematic and cannot be trusted when one deals with the problems where slow and 
fast time scales cannot be separated. This point has been often overlooked, and kinetic 
relations implying a universal elimination of microscopic variables were used in situations 
where both G and V are changing fast as, for instance, in the case of pinning- depinning 
phenomena. 

In what follows we present a detailed adaptation of the above ideas to the case of a 
martensitic phase transition. Martensitic phase boundaries are particularly convenient 
for the demonstration of the main principles of our approach because these plane defects 
may be adequately represented already in one-dimensional models. To emphasize ideas 
we consider the simplest case of a phase boundary with overdamped dynamics. At the 
microscale, the analysis of the non-steady evolution of the core region of such phase 
boundaries requires a study of a dynamical system with an infinite number of degrees of 
freedom. At the macroscale the phase transition is modeled as a singular surface (jump 
discontinuity) whose evolution is governed by a kinetic relation fSB]. 

We pose the question of whether an intermediate description is possible when the inter- 
face is equipped with a small number of "mesoscopic" degrees of freedom whose dynamics 





$(z/; G) = Gp cos u - Gu. 
Then a direct computation gives (e.g. [T^ [T]) 



(3) 




(4) 
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reproduces the main transient effects. An early example of such a reduced description 
can be found in the theory of dislocations where static defects are often represented as 
effective particles in the Peierls-Nabarro (PN) landscape The static FN landscape 
is obtained by relaxing all microscopic variables other than one collective variable inter- 
preted as the macroscopically observable location of the core. The idea of a tilted PN 
landscape has been heuristically apphed to the description of dynamic dislocations in close 
to continuum limit (121 |16l |l7l [IHl [191 ES, EHl US] . This approach, however, cannot be used 
in principle in the strongly discrete case when the dislocation core is atomically narrow 
[T3t [20] . In order to deal with this limit it is natural to abandon the idea of macroscopic 
collective variables and trace instead the dynamics of the particular discrete elements 
while again enslaving all others. Such elements were called "localized normal modes" in 
[32] and "active points" in [23j. In some special cases it has been proven rigorously that 
dynamics of "active points" corresponds to the motion along the center manifold of the 
infinite-dimensional dynamical system |9|. 



In this paper we extend the approach of "active points" from the immediate vicinity 
of the depinning point where it has been formally justified to a broader class of nonsteady 
motions of defects. The key to our method is the assumption that the dynamics of only 
a few bonds located in the core region has to be resolved fully. The other bonds remain 
confined to near bottoms of their respective potential wells, and their small adjustment 
to the changing conditions can be treated as instantaneous. 

The main question confronted by such theory concerns the minimal number of internal 
variables which is sufficient to capture the response of a defect to a particular class of 
external perturbations. In this paper we propose to base the selection of the minimal set 
of internal variables on the careful analysis of the relative "activity" of the variables in the 
traveling wave solution. As a result, we obtain kinetic equations which describe faithfully 
only small deviations from such steady states. Strong deviations from the traveling wave 
ansatz have been rigorously studied for several classes of nonlinear equations in the case 
of periodic obstacles (pulsating traveling fronts, e.g. [TU]). Our approach deals with 
the singular limit of such problems (localization of the fronts), which, as far as we know, 
has not yet been studied in the mathematical literature. 

In our construction we use essentially the fact that the exact traveling wave solution 
of our discrete problem is known in the case of a piecewise quadratic interaction potential 
[lOl [H] . Using this solution as a benchmark, we construct a simple two-dimensional 
dynamical system which in the case of a constant macroscopic driving force generates a 
remarkably good approximation for the kinetic relation originating from the study of the 
full infinite-dimensional dynamical system. We then apply this low-parametric description 
to the case when the driving force is a given function of time and show how it can be 
used to obtain a rheological model of the transformation kinetics which is nonlocal in 
time. We also provide an example showing that depending on the rate of external driving 
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the responses of the phase boundary based on the kinetic relation and on the associated 
kinetic equations can be markedly different. A systematic application of the obtained 
kinetic equations to solutions of the particular macroscopic boundary value problems 
involving strongly inhomogeneous media will be given elsewhere. 

The paper is organized as follows. In Section 2 we formulate the singular macroscopic 
problem which requires a microscopic closure. In search for such a closure we turn in 
Section 3 to a specific microscopic model. To illustrate the existing methods of model 
reduction in relation to moving defects, we construct in Section 4 the static PN landscape 
and then use its tilted version in Section 5 to describe the steady state dynamics of 
phase boundaries. Since this approach is not satisfactory, we develop in Section 6 a 
more systematic one-parametric kinetic equation and introduce a concept of dynamic PN 
landscape. We then find that the one-parametric model is only adequate in the immediate 
vicinity of the depinning point and proceed in Section 7 with the systematic derivation of 
the two-dimensional system of kinetic equations. We find that such description works well 
in the whole interval of admissible velocities. In Section 8 we choose a specific loading 
program and compare the predictions of the nonlocal model with the predictions based 
on the classical kinetic relations. Our concluding remarks are collected in Section 9. To 
make the exposition self-contained, we present in the Appendix a concise derivation of 
the exact kinetic relation. 

2 Macroscopic model 

Let u{x, t) be the one-dimensional continuum displacement field. The macroscopic energy 
of a bar undergoing martensitic phase transition can be written as 



Here Ut = du/dt is the macroscopic velocity, Ux = du/dx is the macroscopic strain, p > 
is the reference mass density, and the energy density (piux) is represented by a double- well 
potential which will be specified later. 

The dynamic equation corresponding to the energy ([5]) is 



where a{ux) = 'P'iux) is the stress-strain relation. On discontinuities, the equation ([6]) 
must be supplemented with the jump conditions. Let /_ and /+ denote the values of /(x) 
to the left and to the right of the interface, and introduce the notations |/] =/+ — /_ for 
the jump and {/} = (/+ + /-)/2 for the average of / across the discontinuity. Then the 




(5) 




(6) 
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parameters on a discontinuity must satisfy the classical Rankine-Hugoniot jump conditions 



{utj + Vlu,,j = 0, pVlutj + = 0. 



(7) 



The entropy inequality, which must also hold on the jump, can be written in the form 



is the configurational (driving) force. In one-dimensional problems imposing the entropy 
inequality can remedy local nonuniqueness only in the case of supersonic discontinuities 
(shock waves). For subsonic phase boundaries one must specify the rate of entropy pro- 
duction as well. This is usually done phenomenologically through the kinetic relation 
G = G{V) between the driving force and the velocity [31 [36l ET]. If the microscopic 
model is known one can instead compute the function G{V) directly from the microscopic 
traveling wave solution (e.g. [33l [311 [Ml SQl E]). 

Behind the assumption that there exists an algebraic relation between the driving force 
G and the velocity V lies the idea that the width of the transition layer is so small that the 
internal relaxation to the traveling wave profile takes place much faster that the variation 
of the driving force. However, the convergence to the sharp interface limit is expected 
to be nonuniform with respect to the parameter describing the rate of configurational 
loading. For instance, if the phase boundary interacts with an obstacle, the structure of 
this obstacle should be sufficiently diffuse for the algebraic kinetic relation to hold within 
certain error bars. If the profile of the obstacle is sufficiently sharp, the time scales of 
internal relaxation and of external driving are comparable and the traveling wave profile 
is not an adequate ansatz. In this case one can expect a uniform double-parametric 
asymptotics in the spirit of [6] to be more complex than is assumed by the computations 
based on the traveling waves. 

In what follows we do not attempt a rigorous asymptotic analysis in the general 
situation and instead choose a simple microscopic model and formally derive a set of 
differential kinetic equations describing the transient response for the inner structure of 
a phase boundary. 

3 Microscopic model 

Consider an infinite chain of particles, connected to their nearest neighbors (NN) through 
viscoelastic springs and to their next-to- nearest neighbors (NNN) by elastic springs (see 



TZ = GV>0 




where TZ is the rate of energy dissipation and 




(9) 
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bi-stable viscoelastic unit 




linear spring 
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Figure 1: The discrete microstructure with viscoelastic nearest and elastic next-to-nearest- 
neighbor interactions. 



Figure [T]). Suppose that in the undeformed configuration the NN and NNN springs have 
lengths e and 2e, respectively. Let Un{t) denote the displacement of nth particle at time 
t with respect to the reference configuration. We associate with the deformation of nth 
NN spring a discrete measure of strain 



Ur. 



Un-1 



For the viscoelastic NN springs we assume the following constitutive relation: 

/nn(w, w) = 0nn(^) + 



(10) 



where ^ > is the viscosity coefficient. To describe domain boundaries the function 
0nn('U^) must be at least a double- well potential; to obtain explicit solutions, we assume 
that this function is biquadratic: 



0Nn(w) 



1 2 



|/i(w — a)^ + ^a[wc 



w < Wc 

W > Wr- 



(12) 



Under these assumptions the NN elastic units can be in two different phases, depending on 
whether the strain is below (phase I) or above (phase II) the critical value Wc- The elastic 
modulus in each phase is > 0, and the parameter a > measures the transformation 
strain. To simplify calculations, we assume that the NNN interactions are linearly elastic: 



fNNN{w) = 2'JW. 



(13) 



Here we defined Wn = {wn+i + Wn)/'2 as the strain in the NNN spring connecting (n-|- l)th 
and {n — l)th particles. 

The dynamics of the chain is governed by the following system of ordinary differential 
equations: 



peiln =/i[ty„+i - Wn - d{wn+i - Wc)a + 6{wn - Wc)a\ 

+ l{Wn+2 + Wn+1 - Wn - Wn-l) + Hw^+l - Wr. 



(14) 



7 



Here p > is the mass density of the chain and 6{x) is the unit step function. To ensure 
stabihty of the chain we require that 

E = fx + A-f>0, (15) 

where E is the homogenized macroscopic elastic modulus; following [381 EH]; we also 
assume that the NNN interactions are of ferromagnetic type, meaning that 7 < 0. In the 
limit e/L —y 0, where L is a macroscopic length scale, we recover ([6]) with the homogenized 
stress-strain law [H] 

cr{w) = E{w — A9{w — Wc))- 
Here A is the macroscopic transformation strain: 

A = f. (16) 

There are two time scales associated with this problem: the time scale of inertia. Tin = 
£:a/ p/E, and the viscosity time scale, T^isc = i/E. In this paper we consider the over- 
damped limit when T^\sc » ^in, i-e. 



e » e^/pE. (17) 
We can nondimensionalize the problem using T^isc as the time scale and letting 

- tE _ Un _ Wn _ Wc 

t = —r, Un = Wn = Wc = —r- (18) 

Dropping the bars on the new variables, we obtain dimensionless system equations: 

Wn - Wn+1 = Cr{Wn+l) - ^{Wn) + D{Wn+2 + ^n+l - Wn - Wn-l), (19) 

where 

a{w) = w-e{w- Wc) (20) 
is the rescaled macroscopic stress-strain law. The dimensionless parameter 

D = -—>0 (21) 

measures the relative strength of NN and NNN interactions. In what follows it will also 
be interpreted as a measure of coupling of the bistable units. 
Observe that system flTI?]) can be "integrated" , yielding 

Wn = D{Wn+l - 2Wn + Wn-l) - ^{Wn) + (T. (22) 
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Here a = a{t) is the time-dependent applied stress. One can see that (1221) also governs the 
overdamped Frenkel-Kontorova model P [HI ESI Elj ; a subtle but important distinction 
of the present setting is that here the coupling coefficient D is independent of e. 

In the overdamped limit the macroscale governing equations reduce to {a{ux))x = 
0, or a{ux) = while the jump conditions take the form = [?T]. The driving 

force (111) is then a function of time: 

G(t) = ait) - au, 

where 

aM=w,~l/2 (23) 

is the Maxwell stress. The function G{t) also enters the microscopic equation (l22ll which 
can be rewritten in the gradient-flow form 

w = -VW(w;G(t)). (24) 

Here w G M°° is the vector of strains, the gradient is taken with respect to w, and 

W = (^wl- {Wn-Wc)0{Wn-Wc) + D{Wn+l-Wnf - {crM + G{t))Wnj, (25) 

n=—oo ^ ' 

is the dimensionless energy of the system. 

Our main goal will be to approximate the infinite-dimensional dynamical system ( IMl) 
by a finite-dimensional reduced dynamical system of the type 

= -a.V(!?{i^;G{t)), (26) 

where ct is the effective mobility matrix. The gradient is taken with respect to the order 
parameter u e M.^ where the integer-valued parameter K defines the dimensionality 
of the reduced system. After the solution of the vector equation fl26l) is known, the 
approximation of the discrete field should be recoverable from the auxiliary relations 
Wn{t) = w„(z/i(t), . . . , z//^(t)) describing the recovery of relaxed "non-order-parameter" 
variables. 

In the absence of a rigorous method condensing to (!26|) we base our formal re- 
duction on the detailed study of a traveling wave solution of (1241) in the form w„(t) = 
w{n — Vt). This solution, which is known explicitly (see Appendix), corresponds to the 
case when G(t) = const and describes a steadily propagating phase boundary. We empha- 
size that in our method the G = const solution is used only to suggest the dimensionality 
of the reduced system, while the actual behavior of the resulting dynamical system is fully 
driven by the time- dependent G{t). 
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4 Classical PN landscape 



The simplest one-parametric {K = 1) reduction of the infinite-dimensional system ([2i 
is based on the idea of Peierls-Nabarro (PN) energy landscape which was first developed 
for dislocations (see the review [29]) and then generalized to phase boundaries in [38] . 
The construction is fully static and therefore restricted to the case of sufficiently small G. 
Since the ideas behind the construction of the PN landscape are used later to obtain the 
approximate kinetic equations, we outline the main steps below; see [38] for more details. 

Assume that G = a — ctm = const and compute stable monotone solutions of the 
equilibrium difference equations fl24l) . For the phase boundary located at n = m we 
obtain (see also [HI [IE]): 

^^^_ exp(A(n-m-l/2)) 

"V ; ^ ^ ^ , exp(-A(n - m - 1/2)) ^ ^ 

2cosh(A/2) ' 

where 

X{D) = arccosh (^1 + . (28) 

The admissibility constraints 

> Wc for n < m, w'!^ <Wc for n > m + 1 (29) 

determine the constraints on G; the set of driving forces satisfying these constraints 

constitutes the trapping region. One can show that for D > Q the strain profile ( 1271) is 

monotone, so the constraints (l29l) can be replaced by > Wc and w'^j^^ < Wc- The 
trapping region is then given by 

\G\ < Gp, (30) 

where 

Gp(D) = - tanh - = , ^ (31) 

is the desired expression for the Peierls threshold [HI |38] . 

One can see that at IGI < Gp the phase boundary can be in an infinite number of stable 
equilibrium configurations fl27|) parameterized by m. In order to connect these stable 
states one needs to consider a path involving non-equilibrium intermediate configurations. 
For instance, take two equilibrium configurations, one with the phase boundary located 
at n = i — 1 and the other one at n = i Suppose that phase II is located behind the phase 
boundary and consider all paths connecting the two equilibria along which the ith NN 
spring is the only one that changes phase, while all other springs stay in their respective 
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energy wells. As the ith NN spring crosses the critical value of strain Wc the system goes 
through an energy barrier and the next task is to select the path that involves the minimal 
barrier. 

To this end we fix Wi and minimize the energy of the chain with respect to all strains 
Wk with k ^ i. We obtain: 

{l + 2D)wk-D{wk+i + Wk-i) =crM + > i + 1 ^^^^ 



Solving equations (l32l) at k < i — 1 and k > i + 1 and requiring that the corresponding 
solutions, when extended to k = i, both equal Wi, we obtain the following strains along 
the path connecting the two equilibrium statesQ 

^ / G + l + e^(^-)K-aM-G'-l), k<i 
w,-aM + ^^^ e-^(^-)K -au-G), k>i. 

At the saddle point Wi = Wc the energy reaches its maximal value; due to the nonsmooth- 
ness of the biquadratic potential at w = Wc, the corresponding states are singular. 

The next step is to construct a global path that connects not just two but all equivalent 
equilibrium points. To this end we replace the order parameters Wi which were operative 
in the consecutive segments of the path by a global order parameter u defined implicitly 
by 

=w*{u) = aM + G + ^^^^ + " M) tanh(A/2). (34) 

Here [i'] denotes the integer part of u. At the integer values of u we have W[^] = wl, and 
as u approaches [u] + 1 from below, W[^j linearly increases to Wu. The function w*(z/) is 
periodic with period 1 and has jump discontinuities at v = One can see that v can 
be interpreted as the macroscopic location of the interface. 

Using and recalling fl33l) . we obtain strains along the path that are now parametrized 
by v: 

G+1 + e^('=-M)(w;*(z/) - ctm - G - 1), k<[u] 



Wk-crM+^^_^ ^_A(fc-M)(^*(^) - aM - G), k > [u]. ^^^^ 

To evaluate the energy W along the path fl5^ we first renormalized it as ^{u-jG) = 
W{iy) - W(0). Substituting ([35D in ([25]), we obtain (recall that < G < Gp) : 

$(//; G) = Gp{u - H)2 - G[u] - (G - Gp + 2{u - [u])Gp)9(u - [u] - ■ (36) 



^Notice that along the selected path the order parameter Wi increases from its value in the first 
minimum, wl = u'*^^(G') = (TM + G + exp(— A/2)/(2cosh(A/2)) < Wc, to its value in the second minimum, 
wu =wl{G) = (TM + G + f -exp(-A/2)/(2cosh(A/2)) > Wa- 



ll 



G = 1.2Gp,K=l 



Figure 2: One-dimensional energy landscape $(z^;G') at D = 0.5 and different G. The 
classical static PN landscape corresponds to G < Gp. The dynamic potential above the 
Periels threshold G > Gp was constructed using the K = 1 approximation. 



Observe that the minima of the resulting PN landscape $(z/; G) are located at the integer 
values of u and correspond to stable equilibrium states, while the singularities at u = Ui 
represent unstable equilibria (saddle points) where Wiivi) = Wc (see Fig. [2]). The energy 
barriers separating equilibria at = i and i> = i + 1 are independent of i: 

$(..) - $(z) = - 



4Gp 

At G = Gp the stable equilibria become marginally stable; see the curve G = Gp in 
Fig. [21 Above the Peierls threshold (Gp < G < Gs), the stable equilibria cease to exist, 
and the above construction of the PN landscape becomes invalid. 



5 Dynamics on a tilted PN landscape 

One very heuristic but rather common way to extend the idea of the PN landscape to 
dynamics is to take the classical PN landscape at G = and tilt it to account for nonzero 
G. Then one can obtain an approximate kinetic relation by considering the overdamped 
motion of a material point on this tilted landscape and studying the dependence of the 
average velocity on the driving force [7]. 

The tilted PN landscape can be defined as follows 

<I>t(z/;G) = $(// + 1/2; 0) -Gz/- (37) 

where we have added a constant to ensure that $t(0; G) = 0. It is not hard to see that 
at G < Gp the energy (!36|) agrees with (1371) up to the appropriate shifts in horizontal and 
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(a) 




tilted, 
G=1.2Gp 




Figure 3: (a) The tilted energy landscape (grey) and the DPN landscape for K = 1 (black) 
at D = 0.5. The two landscapes coincide at G = Gp. (b) The exact kinetic relation G{V) 
from the Appendix (solid curve) and its approximations via the tilted landscape (grey 
curve), K = 1 approximation (dotted curve) and K = 2 approximation (dashed curve). 



vertical directions: 



<^iu+ g] - : = $t(^; G), 0<G<Gp. 



2Gf 



4Gp 



Note that the two landscapes coincide exactly only at the Peierls threshold. Indeed, in 
the actual PN landscape (|36l) the local energy minima are always located at the integer 
values of u, reflecting the fact that in a discrete system position of a phase boundary always 
coincides with a lattice point. As G approaches the value Gp the points of local maxima 
of the PN landscape also move towards integer values and eventually the minumum and 
maximum points merge. Meanwhile, in the case of the tilted PN landscape fl37|) the local 
maxima are always fixed at the integer values, and the minimum points move instead. 

Despite these deficiencies, the tilted PN landscape has a strong advantage: it can be 
formally extended beyond the Peierls threshold G = Gp. One can then try to model 
the dynamics of a phase boundary by the overdamped motion of a particle in the tilted 
landscape (see the grey curve in Figure [3^1). The equation of motion takes the form 

z> = -a<^'^{u; G) = -a(2Gp(z/ - [u] - 1/2) - G). 

Solving this equation subject to the condition z/([Vt]/V) = [Vt] yields 

^(t) = [Vt] + ( 1 - e-2"Gp(*-[v*]/v^)^ . 

2Gp 
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Imposing the periodicity condition z/(([\^t] + l)/V) = [Vt] + 1, we obtain the desired 
kinetic relation 

1 I ~2aGp/V 

GiV) = gP i_,-2.Cp/v- (38) 

As expected, it gives the correct hmit G Gp when V 0. However, in contrast 
to the exact kinetic relation obtained in the Appendix (see (185|) ) the function (!38|) is 
unbounded: at large V it grows as GiV) ^ V/a (see Figure [3b)- In addition, since the 
physical meaning of the dynamic variable i/(t) is obscure, it is not clear how one can 
recover the strains Wn{t). We therefore abandon this heuristic path and search for a more 
adequate representation of the actual dynamics of the chain. 



6 Reduced model with one internal variable 

We now take a more systematic point of view, choose a single variable and follow its 
exact dynamics while minimizing the energy with respect to all remaining variables. The 
resulting model is similar to the "single-active-site theory" of [9| [23]. 



6.1 K = 1 model 

Suppose that G = const and assume that the motion in the reduced system is periodic. 
Assume also that during each period dynamics of only a single NN spring, located right 
behind the phase transition front and actually changing the energy well during this period, 
needs to be traced in all the detail, while all other springs can be "enslaved". We can 
always assume that at t = the 0th NN spring has just switched to phase II: Wo(0) = Wc- 
Therefore, during the time period < t < 1/V the "active strain" is Wq. Minimizing the 
energy with respect to all other strain variables, we obtain 

I [Wq — Wc — G + l/2)e n > 1. 



Substituting these expressions in (1221) for n = 0, we obtain a single equation governing 
the dynamics of the active point: 

Wq = -wo{l + 2D{1 - e-^)} + Wc + G+l/2 + 2D{G + Wc){l - e'^). (40) 



6.2 K = 1 kinetic relation 

To see whether the ensuing dynamics is compatible with the exact kinetic relation 0851) 
in the Appendix, we solve (HUl) subject to the boundary condition Wo{0) = Wc, obtaining 

Wo{t) =Wc + {G + Gp){l- exp(--^] ], (41) 
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Then, setting = 1 in (]39l) . we find 

wi{t) =Wc + G-Gp- e'^{G + Gp) exp 



2Gv 



The second boundary condition wi{l/V) = Wc yields the desired approximation of the 
kinetic relation: 

^^p(2Gpk) ~ exp(-A) 

Clearly, G{V) — > Gp as — 0, and in agreement with ( l85l) G'(V^) tends to the spinodal 
value Gs = 1/2 as goes to infinity. Note, however, that the asymptotic behavior at 
small V, G — Gp ~ exp(— g^^), differs from the exact asymptotics 



obtained in ITT]. 



6.3 K = 1 dynamic PN landscape 

To construct the dynamic PN landscape we need to to patch different segments of the 
dynamic trajectory into one by introducing a single order parameter. Recall that during 
each time period k/V < t < {k + ^)/V the active strain is Wk{t), and the corresponding 
solution can be obtained by replacing n in the expressions above for k = hj n — k and 
t hj t — k/V. Introduce a new continuous monotonically increasing function 

1 - exp(-5J^(« - 1^)) 

i expi^ 2Gp>' (44) 
G + Gp f I f [VtV ' ~ 



(Gs-G)(e^-l)i "V 2GpV V 

where [z/] = [Vt] and note that z/(t) takes consecutive integer values at the beginning 
of each time period. Therefore, it is exactly the dynamic extension of the parameter u 
which we used to obtain the static PN landscape. In terms of this new "global" order 
parameter, we obtain the following expression for an active strain W[u]{t) (over the time 
period [u]/V < t < {[u] + 1)/V): 

wi.]it) =w, + (e" - l)(Gs - G)(z/(t) - [z/(t)]). 
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(a) 



(b) 





Figure 4: (a) Kinetic relations G{V) for infinite-dimensional dynamics ( l85i) (solid curve), 
its small- velocity approximation f H3l) (grey curve) and K = 1 approximation (dotted 
curve), (b) Comparison of strain trajectories obtained in the K = 1 reduced model 
(dotted curves) and in the original infinite-dimensional model (solid curves). Parameters: 



V = 0.1, D = 0.5, Wc = 1. 



Time evolution of all strains can be then written in terms of u a^ 

Wnit) =Wc- 1/2 + G+ 

1 + ((e^ - l)(Gs - G){u{t) - [Vt]) -G- Gs^e^^""!^*!), n < [Vt] (45) 
(e^ - l)(Gs - G){u{t) - [Vt]) -G + Gs^e-^("-[^*]) n > [Vt] + 1. 



Now we are in a position to construct the global DPN potential ^{u; G) which serves as 
a landscape for the overdamped dynamics of t^it): 



-«$'(//; G). 



(46) 



The potential $(i/; G) is piecewise quadratic and we can select the coefficient a = l/(4Gp) 
to ensure that the function (HH) is the solution of (HSl). As the result, we obtain 



$(//; G) = Gp 



V 



2(G + Gp) 



(eA_i)(Gs-G) 



(47) 



^Although v{f) is by construction a continuous function of time, the strain variables have jump dis- 
continuities at t = njV ^ n = 1,2 ... . This is due to the fact that our approximation treats non-active 
strains as if they were equilibrated and neglects their real dynamics. As a result, at the end of each time 
period when we switch to the new active strain, the old one has to increase its value discontinuously in 
order to reach equilibrium with the new active point. The associated jumps can be computed explicitly: 
Klt=(„+i)/y = K]t=i/y =w'o(l/'^" + 0)-u;o(l/'^"-0) = 2sinhA(G-Gp). 
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Notice that this potential coincides with the equilibrium PN landscape at G = Gp (see 
Figure ED. 

In Figure |Dd we compare the one-dimensional dynamics fl46p with the full infinite- 
dimensional dynamics at G = const (see Appendix). As expected, the evolution of the 
active strain is captured quite well, but there is a visible deviation from the actual values 
for the strains whose adjustment was assumed to be instantaneous. The corresponding 
kinetic relations are compared in Figure UK- One can see that the reduced model with 
K = 1 provides a good quantitative approximation of the exact kinetic relation until 
V = 0.2. At higher velocities it deviates substantially from the exact relation, although 
both tend to the same spinodal limit Gs at infinite V. The observed discrepancies in the 
case G = const suggest that the reduced model with i^' = 1 is oversimplified and cannot 
be used as the base for the construction of kinetic equations. 



7 Reduced model with two internal variables 

To obtain a reduced model which works well in the whole range of admissible velocities, 
we need to consider more "active points" in the core region of the defect. From Figure Hb 
one can see that in each period the natural choice for the expanded set of order parameters 
would be Wq and Wi. 

7.1 K = 2 model 

Assume again that G = const and minimize the energy with respect to all strain variables 
other than Wq and We obtain the following recovery relations 

( We + G + l/2 + (w;_i-w, -G-l/2)e^("+^), n < -1 
Wn= I Wo, n = (48) 

[ Wc + G - 1/2 + (wi - Wc - G + l/2)e^(i~") n > 1. 

Observe that the dynamics of the variables W-i{t), wo{t) and wi{t) in ( l48l) is not inde- 
pendent. Indeed, the three active points satisfy the dynamic equations 

w^i = D{wq — 2w_i + w_2) — w^i + Wc + G + 1/2 
Wo = D{wi - 2wo + w-i) - wq + Wc + G+1/2 (49) 
Wi = D{w2 — 2wi + Wq) — Wi + Wc + G — 1/2 

•^Notice that instead of minimizing out the "enslaved" variables at each value of G we could have 
also used the tails of the known traveling wave solution (see Appendix). However this leads to very 
complicated implicit algebraic relations and in the interest of transparency we decided not to pursue this 
more rigorous approach in this paper. 
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Substituting the expressions for W2 and W-2 from (l48i) in ( H9i) . we obtain 

= {-2D - 1 + e-^D)w_i + Dw^ + {wc + G + 1/2)(1 + D{1 - e"^)) 
?i;o = {-2D - 1)wq + Dwi + Dw^i + Wc + G + 1/2 (50) 
= {-2D - 1 + e~^D)w^i + Dwq + (w^ + G - 1/2)(1 + D{1 - e"^)) 

If we now compare the equations governing the dynamics of W-i{t) and Wi{t), we see that 
they differ only by a constant term in the right hand side. This allows us to reduce (l50i) 
to a two-dimensional system for 

/X /X /X W_l(t)+Wl(t) 

x(t) = u;o(t), y(t) = ^^i^-^^. (51) 

Note that 

^TiW = ± 2(l + D(2-e-A)) = ± ' 

where we used ( l28i) to obtain the second equality. The two variables: x(t), describing the 
dynamics of the transforming spring, and y{t), describing the average strain in the core 
region, must satisfy the following system of equations: 

x = {-2D-l)x + 2Dy + Wc + G + l/2 (52) 

y = Dx + {-2D - 1 + e-^D)y + {w^ + G')(l + D{1 - e"^)). (53) 
This system will be used as an approximation of the full infinite-dimensional dynamics. 

7.2 K = 2 kinetic relation 

We first check that the traveling wave solutions of the original system (see Appendix) can 
be reproduced by the reduced system fl52l) . fl53|) . Such solution must be subjected to the 
following constraints. First, we must require that 

x{0) = w,, (54) 

a condition that ensures that the 0th NN spring (recall that x = Wq) has just transformed 
to the new phase at t = 0. Second, we require that 

Wi{l/V) = Wc, 

or ^ 

y{l/V)=W,+ ^—^. (55) 

18 



(a) (b) 




Figure 5: (a) Kinetic relations G{V): infinite-dimensional system (185|) (solid curve), its 
small- velocity approximation (H3l) (grey curve) and K = 2 approximation (dashed curve); 
(b) Comparison of the strain trajectories obtained in the K = 2 model (dashed curves) 
and in the infinite-dimensional model (solid curves). Parameters: V = 0.1, D = 0.5, 

Wc = 1. 



This means that ai t = 1/V the first NN spring reaches the critical strain, marking the 
end of the period. Finally, periodicity requires that 

Using ( l48l) at n = 2 together with ( l55l) . we can see that this boundary condition reduces 
to 

wiiO)=w,~iGs-G)il-e-^), 

or 

yiO) = w, + 0(1-6-^). (56) 

Next, by solving the system (l52l) . (153!) subject to the initial conditions (l54l) and (!56l) . we 
obtain 

x{t) =Wc + G + Gp + cie'i* + C2e"2* 

, ; (57) 

y{t) = Wc + G + e-^Gp + — (ci(l + Vl + 8e^)e''' + C2(l - Vl + 8e2^)e''^*), 



where ^ 

^i-2 = -7—: ^^7T7^(l + ^e-2\lT VTTS^)^ (5^ 



4sinh"^A/2) 
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Figure 6: The error in the approximation of the the exact kinetic relation Gex(^) by the 
reduced models with K = 1 and K = 2 generating approximate kinetic relation Gap(^)- 



(note that r2 < ri < 0) and 

(3 + Vl + 8e2^)(2G' - 1 + e^(l + 2^)) 

4(l + e^)Vl + 8e2^ 
(Vl + 8e2^ - 3)(2G - 1 + e^(l + 2^)) 



Cl 



C2 



4(1 + e^)Vl + 8e2^ 



The application of the boundary condition ( l55l) yields the desired approximation of the 
exact kinetic relation: 

Gp + ({1 + 2e2^ + Vl + 8e2^)e'^i/^ + (Vl + 86^^ - 26^^ - l)e^2/v^ 

G'(y) = 4(l+e^)V l+8e2^ ^ / 



2A/l+8e 



(^(1 + 2e2^ + Vl + 8e2^)e^'i/^ + ( VTTS^ - 2e^^ - 



(59) 

First we notice that the approximate kinetic relation (!59|) satisfies the constraint 
G(0) = Gp. Then, as V tends to infinity, G{y) ^ Go, = 1/2 with the asymptotics 
G — Gs ~ exp(— |ri|/V^). The global comparison of the approximate and exact (see Ap- 
pendix) kinetic relations is presented in Figure [51 One can see that in view of how few 
degrees of freedom are involved in the approximation, the agreement is remarkable, at 
least for small (strongly discrete limit) to moderate values of D (see Figure [6]) 



*At larger D (strong coupling limit) the core region delocalizes and more and more active points need 
to be included [9] . Another possibility is to change the type of the approximation from "active points" 
to "collective variables" , which is more in tune with a close to continuum character of the model in this 
limit. We will not pursue these ideas in the present paper. 
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(a) (b) 




Vl 



Figure 7: (a) Two-dimensional energy landscape $(z/i,z/2). (b) Level sets of the energy 
along with the path of the effective particle i^{t). Parameters: V = 0.1, D = 0.5, Wc = 1. 

7.3 K = 2 dynamic PN landscape 

We can now introduce the global order parameter and reformulate the G = const dynamics 
in the form of a gradient flow on a fixed two-dimensional DPN landscape. Recall that 
expressions f l571) . flSTj) and fHHj) define the strain trajectories w„(t) only over the time 
period < t < 1/V. The extension of this solution for k/V < t < {k + 1)/^ is obtained 
by replacing n in the above formulas by n — k and t by t — k/V. To patch together different 
periods we need to introduce a set of two order parameters that change continuously and 
monotonically with t. We denote these variables as t'i(t) and i'2{t) and define them 
by the conditions [t'i(t)] = ['^2(^)] = l^t]. Then the variables x(t) and y(t) can be 
extended periodically to any t > as follows: x(t) = Wc + {i^i(t) — [Vt]){x{l/V) —Wc) and 
y{t) = y{0) + (z>2(t) - [Vt]){y{l/V)-y{0)). Under the assumption that V = T{G) , where 
T{G) is the inverse of GiV) from (15^ . the dynamics of the vector field i>(t) = (t'i(t), i^2it)) 
can be represented in the form 

= -ctV^{i>;G). (60) 

Here we introduced a mobility matrix a and a two-dimensional effective DPN potential 
$(z>,G). 

We can further simplify the system ( l60i) by diagonalizing the mobility matrix a. 
Observe that in fl58|) we have \r2\ > \ri\ > 0, so that the eigenvector (^(A/l + Se^+l), 1) 
corresponding the eigenvalue r2 is the slow direction, while the eigenvector (^^-(l — 
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VT+~8e^), 1) that corresponds to ri is the fast direction. Introduce the slow variable 

l-exp{n{t-[Vt]/V)) 



and the fast variable 



Mt) = [vt] 



Mt) = [Vt] + 



1 — exp(ri/l^) 
l-exp{r2{t-[Vt]/V)) 



(61) 



(62) 



1 - exp(r2/V') 

In terms of the new variables, the time evolution of all strains is given by 

w, + G + 1/2 + {Gp - 1/2 + ici(l + Vl + 8e2^)(l - (1 - e-'/^)iiy,{t) - [Vt])) 

+C2(1 - Vl + 8e2A)(l - (1 - e'^/^){iy2{t) - [W]))}e^("-[^*1), n < [Vt] - 1 

w, + G + Gp + ci(l - (1 - e'-i/^)(z/i(t) - [Vt])) 

+C2(1 - (1 - e'-^/^)(z/2(t) - [Vt])), n = [Vt] 

1 



Wn{t) 



w, + G -1/2 + {Gp + 1/2 + -ci(l + Vl + 8e2A)(l - (1 - e'-i/^)(z/i(t) - [Vt])) 
+C2(1 - Vl + 8e2^)(l - (1 - e'-^l^){u2{t) - [W]))}e-^("-[^*]), n > [Vt] + 1. 



The vector field ^{t) = (i^i(t), i^2(^)) satisfies 

iy = -aV<I>(i/;G), 



(63) 



(64) 



where the mobility matrix a = diag(ai, ^2) is now diagonal (here ai > and 02 > are 
constants). The potential $(z/i, 1/2; G) is piecewise quadratic; it depends on V and hence, 
through the kinetic relation V = T{G), on the driving force G. Using fl6T|) and fl62l) . we 
obtain 



$(z/i,z/2;G') 



ri 



2ai 
'202 



>2 - h])^ 



exp(ri/V(G)) 
2 

exp(r2/V(G)) 



T^i + [J^i] 
1^2 + [1^2] 



To determine the effective viscosities a\ and 02, we require that at Gp the DPN 
potential equals the relative Gibbs free energy of the system: 



$(iy; Gp) = W{u] Gp) - >V(0; Gp). 



This yields 



(Vl + 8e2A + 3)Gp' 



a2 



riVl + 8e2^ 



(Vl + 8e2^ - 3)Gp 



(65) 



22 



and, finally, 



(Vl + 8e2^ - 3)Gp r, ^ 2 . 



(66) 



The dynamic PN landscape (DPN) $(z/i, z/2; G), along with the trajectory of the effective 
particle, is shown in Figure [3 The comparison of strain trajectories in the two-dimensional 
reduced theory with the exact result at G = const (see Appendix) shows the considerable 
improvement over the one- dimensional approximation. For instance. Figure [5b compares 
the evolution of strains near the phase boundary over the first three time periods at 
V = 0.1. While in the case K = 1 the evolution of only the transforming element is 
followed closely over each time period, in the K = 2 case the dynamics of the main 
nontransforming elements forming the core region of the defect is also captured extremely 
well. It is then natural to use the K = 2 approximation as the basis for the construction 
of the nonlocal rheological relation on the phase boundary. 



8 Kinetic equations 

In the previous section we saw that the traveling wave solution of the original infinite- 
dimensional system can be reproduced rather faithfully by the reduced two-dimensional 
system (!52|) and (1531) . This suggests that the K = 2 model may be extended to deal with 
the situation when the driving force is not constant but varies sufficiently slowly. In this 
more general setting we obtain the system 

X ={-2D - l)x + 2Dy + w^ + G{t) + 1/2 

y =Dx + {-2D - 1 + e-^D)y + [w^ + G{t)){l + D(l - e"^)). ^^''^ 

Equations 0671) can be interpreted as the system of differential kinetic equations gener- 
alizing the kinetic relation fl59|) . According to the description of kinetics implied by the 
system fl671) the adjustment of the velocity V{t) to the dynamic configurational loading 
G{t) is not instantaneous and cannot be described by an algebraic relation if the func- 
tion G{t) changes sufficiently fast. Instead, the relationship between the velocity and the 
driving force is expected to be history-dependent: 

V;(t) = J^{G(r),r <t}. (68) 

The goal of this section is to reconstruct the nonlocal relationship (1681) from the system 
(|67|1 and to illustrate the effect of nonlocality. 
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Suppose that the phase boundary is located at n = z at time t = U, meaning that 
in this instance the ith spring has critical strain: Wi{ti) = Wc- As before, we can define 
x(t) = Wi(t) and y{t) = +Wj_i(t))/2 when ti < t < tj+i. To obtain the average 

velocity Vi we need to solve the system of kinetic equations fl67j) in the interval 
under the assumption that G{t) is known. 

Observe first that conditions 

x{ti) = Wc (69) 
and ^ 

y{ti+i) =Wc+ , (70) 

which are the extensions of (IMIl and (!55|) to the interval ensure that the system 

holds from t = ti until the time tj+i when the next spring reaches the critical strain: 

While for given tj+i conditions fl69|) and (I70|) are sufficient to find a unique solution 
of fl67|) . the moment tj+i when the consecutive spring switches from one energy well to 
another remains unknown. If G{t) varies sufficiently slowly, it is reasonable to assume 
that the motion is close to periodic in the sense thacl 

Wiiti) = Wi+i(tj+i). (71) 

This condition allows us to close the system and find tj+i from 

yiU) =Wc + GiU+i)il-e-''). (72) 

To this end, we solve ( 167|) subject to ( l69i) and ( 1721) . obtain for ti <t < tj+i, 

x{t) =g{t) + Ae^^(*-^) /i(r) + e'-^^'-^^ h{T))dr 
Jti 

yit) =h{t) + ^ f {{l + Vl + 8e2A)e^^(*--)A(r) ^^^^ 



+ (1 - Vl + 8e2^)e"^(*-")/2(r)}cir, 



'''This condition is obviously not exact because G{ti) ^ G(ti+i). By using the knowledge of the 
exact structure of the traveling waves at both G{ti) and G(t,;+i) this approximation can be in principle 
improved. However, in the interests of transparency, we do not pursue this more rigorous approach in 
this paper. 
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where we defined 

Vl + Se^^ [2 \ / 



2 



flit) =^V^=^K + - + + + + D(l - e-^)) 



a/1 _|_ 8p2A _|_ 1 1 OpA 

(74) 

We can now determine the strains Wi{t) = x{t) and Wizfii{t) = y{t) ± (1 — e^'^)/2. The 
other strains can be recovered as before, from the constrained energy minimization: 



Wn{t) 



Wc + G{t) + 1/2 + {wi^i -Wc- G{t) - l/2)e^("~^+i), n<i-2 
Wc + Git) - 1/2 + {wi+i -Wc- Git) + l/2)e^(i-'^-^) n>i + 2. 



Condition (!70|) apphed to y(t) in (!73l) yields the nonlinear equation for ti+i. 



—A ,o^A rti+i 



1- e'^ e 



4 



h{U+,) -Wc- ^— — + — / {(1 + Vl + 8e2A)e'^i(*'+i--)/i(r) 



(75) 

+ (1 - a/1 + 8e2^)e'^2(*»+i~")/2(r)}rfr = 0. 



As we have seen in the previous section, at G = const this equation a unique solution 
and one can expect that at least for the case of slowly varying driving force G{t) the 
solution tj+i also exists and is unique. Once it is found, we can obtain the average 
velocity corresponding to the time interval in the form 

V^ = (76) 

Using this iterative process, one can generate the sequence of successive time intervals 
[ti, ij+i] and find the average velocities associated with each of them. This gives a piecewise 
constant function V{ti), which appears at the macroscale as a continuous function V^{t). 
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Figure 8: (a) Time-dependent driving force (1771) . with ^ = 0.1 (dashed curve) and C = 10 
(sohd curve), (b) The corresponding difference between the coordinate s{t) of the phase 
boundary obtained using the instantaneous kinetic response V{t) = r{G{t)) and the 
coordinate s^.(t) calculated using the kinetic equations fl67j) . Here Gq = G{0.3) ~ 0.294, 
r = 0.001 and = 20. 

To illustrate this procedure consider a specific time-dependent driving force of the 
form 

G{t) =Go-S tanh (77) 

Here Gq = G{t^) is the constant background driving force, 6 is the amplitude of the diffuse 
jump in the driving force, and ( is the characteristic time over which the driving force 
decreases by 26. Such perturbation of the driving force may represent an interaction of 
the phase boundary with an extended obstacle. 

We now compare two responses of the steadily moving phase boundary to a perturba- 
tion ( 177|) . One is the instantaneous response governed by the kinetic relation (l59i) . In this 
case we obtain V(t) = T{G(t)) and the location of the phase boundary assuming s(0) = 
can be found by integration from 

s{t)= fr{GiT))dT. 

Jo 

The second one is the history-dependent response governed by the kinetic equations fl67j) . 
In this case the function V^{t) is obtained from the interpolation of the solution of the 
sequence of the discrete problems (175]) . The location of the phase boundary is given by 

i 

3=0 
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The time- dependent separation between the two boundaries both located at the same 
point at time t = but then driven by two different kinetic models is shown in Fig. [HI 
As ( oo, the two kinetic models converge and the difference s{t) — s,,(t) tends to zero. 
However, as Fig. [8] shows, at finite ( the prediction of the model based on the kinetic 
relation may be markedly different from the one based on solving the kinetic equations. 

9 Conclusions 

In the classical macroscopic continuum description of lattice defects the core regions are 
represented by singularities, and the dynamics of these singularities is governed by alge- 
braic kinetic relations. Such description is too coarse to capture the details of the dynamic 
response of the cores of the defects to relatively fast changes of the macroscopic driving 
forces. In particular, it averages out the details of the intricate interaction of the core re- 
gions with localized inhomogeneities which manifests itself through macroscopic velocity 
oscillations and a characteristic acoustic emission. 

To describe the "breathing" of the core regions as well as other transient effects that 
are usually neglected, we proposed in this paper to replace the algebraic relations between 
the macroscopic velocity and the corresponding driving force, which imply internal steady 
state, by a system of differential kinetic equations. The explicit time dependence of 
the driving force makes these equations nonautonomous and allows them to represent a 
dependence between the current value of velocity and the history of the driving force. 
In the constitutive sense the proposed kinetic equations describe a relation between the 
thermodynamic force and the corresponding flux which is nonlocal in time. The ensuing 
rheological description of the moving defect can be qualified as a differential rate model. 

The main idea of our construction is to follow the exact dynamics of only few discrete 
degrees of freedom. The adiabatic elimination of the infinite number of the remaining 
variables outside the core region serves as the matching condition with the classical con- 
tinuum description outside the defect. The choice of the dimensionality of the reduced 
system remains heuristic and is based on the comparison of some special periodic solutions 
of the finite-dimensional system with the discrete traveling wave solutions of the micro- 
scopic infinite-dimensional model. We show that in contrast to the standard approach 
that involves only the center of mass of the defect, the minimally adequate approxima- 
tion includes not only the location but also the internal configuration of the core. Such 
extension of the scope is necessary even if both the macroscopic driving force and the av- 
erage velocity of the defect are constant because the core region may experience periodic 
configurational deformations instead of translating as a rigid body. 

In this paper we presented only the contours of this approach and provided only the 
most elementary illustrations. More systematic study is needed to formulate the formal 
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mathematical structure of our center-manifold-type reduction procedure and to evaluate 
the errors of the approximation in mathematical terms. Other problems include transition 
from overdamped to underdamped dynamics and the generalization for higher dimensions. 
A specific development is also needed to move from a prototypical case of martensitic phase 
transitions to the cases of dislocations, cracks, and more complex singularities. 

Appendix: Traveling wave solution 

To construct a traveling wave solution of the system fl22l) with o = const, we assume that 

Wn{t) = w{r]), r] = n-Vt, (78) 

where V is the dimensionless velocity of the front, which represents a moving phase 
boundary. Since we seek a description of an isolated phase boundary that leaves phase II 
behind, we require that 

w{ri) < Wc for r] > 0, wirj) > Wc for < (79) 

Under these assumptions, fl2^ reduces to 

Vw'{t]) - D{w{r] + 1) - 2w{r]) + w{r] - 1)) + w{7]) = 9{-t]) + a (80) 

At infinity the solution must tend to uniform-strain equilibria of f|T9|) : 

w{ri) — ^ w± as r] ^ ±oo. (81) 

Finally, for consistency, we must also require that 

w{0) = Wc. (82) 

Since the equation (IHOl) is linear, we can be solve it using Fourier transform (see 
[91 [m E] for details). We obtain 



wirf) 



^+1+ E , , forr/<0 



^ ~ Z] , . /, for 7] > 

fcg5+(y)/cAfe(A;,\/) 



where S^iy^ = {k : A{k, V) = 0, Imk ^ 0} are the sets of roots of the dispersion relation 

A(fc, V) = 1 + AD sm^{k/2) - Vik = 0. 
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Figure 9: Kinetic relation at different values of D. 



Continuity of ^(r]) at r/ = gives the relationship between the applied stress and velocity 
of the traveling wave: 



0" = Om 



1 



1 



2 ' ^ kkAk.V) 



E 



1 



""^^ 2 ^ kAk(k,V)' 



Since the difference between applied and Maxwell stresses is equal to the driving force 
G = a — aM (see |lQ]), we obtaiq§ 



G{V) 



1 
2 



E 

fce5+(y) 



A;Afc(A;,V)' 



^5) 



The structure of the kinetic relation fl85l) is illustrated in Figure [9] for different values of 
D. One can see that a.t D = (no NNN interactions) the driving force must be constant 
and equal to the spinodal value G = Gs = 1/2; at = it can take any value between 
and Gs- In the limit D ^ oo the driving force becomes equal to zero at all values of V. At 
large V complex roots in the set S~^{V) tend to infinity and the kinetic curves approach 
the common limit. For D > 1/12 the large- velocity discrete kinetics is well approximated 
by the formula [H] 

V 

G{V) = =, 
2^V^ + 4D -1/3 

which describes the kinetics of the overdamped viscosity-capillarity model governed by 
the partial differential equation Wt = {D — l/12)wj:x — a{w) + a. 



^This semianalytic expression remains implicit mrtil the exact locations of the roots k = k(V) of the 
dispersion equation are known [41 . 
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